Toroidal Oscillations of slowly rotating relativistic star in tensor-vector-scalar theory 
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We examine the toroidal oscillations on the slowly rotating relativistic stars in tensor-vector- 
scalar (TeVeS) theory with the Cowling approximation. As a result, we find that perturbation 
equations describing the toroidal oscillations are same equation form as in general relativity (GR). 
Although the frequencies of toroidal oscillations in TeVeS are not so different from those in GR, the 
momentum inertia depends strongly on the gravitational theory. Thus, observing the frequencies of 
toroidal oscillations and momentum inertia with high accuracy might reveal the gravitational theory 
in the strong- field regime. 
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Q\ 1 I. INTRODUCTION 



rj"| | Via many experiments, the validity of general relativity (GR) has been shown in the weak gravitational field such 
as solar system, while the tests of gravitational theory in the strong-field regime are still very poor. Accordingly, 
the gravitational theory in the strong gravitational field has not been constrained observationally. However, with the 
,£3 • development of technology, it is becoming possible to observe the compact objects with high accuracy and via these 
Qh' observations it could be possible to test the gravitational theory in the strong gravitational field jlj. In practice, some 
possibilities to distinguish the gravitational theory are suggested. For example, using the surface atomic line redshifts 
[2| or gravitational waves radiated from the neutron stars [3j, it could be possible to distinguish the scalar-tensor 
theory proposed in Q from GR. Additionally, Corda suggested the definitive test for GR with gravitational waves 

The tensor- vector-scalar (TeVeS) theory, which is originally proposed by Bekenstein Q, is one of the alternative 
gravitational theory and attracts considerable attention. This theory is covariant formalism for modified Newtonian 
dynamics 0, Hi ■ The reason why this theory gets attention is possible to explain the galaxy rotational curve and 
' Tully-Fisher law without the presence of dark matter Q . Te VeS is also successfully to explain the strong gravitational 
lensing [9( as well as the galaxy distributions through an evolving Universe without cold dark matter [10]. It should 
be noticed that the bullet cluster 1E0657-558 might be good candidate to make a constraint in the gravitational 
theory observationally [ll|, but so far no one has tried this attempt for TeVeS yet. Drawing attention to the strong 
~ . ■ gravitational region of TeVeS, Giannios found the Schwarzschild solution [12 ]. Sagi and Bekenstein found the Reissner- 
Nordstrom solution [13J, and Lasky et al. derived the Tolman-Oppenheimer-Volkoff (TOV) equations in TeVeS and 
produced the static, spherically symmetric stellar models in TeVeS [14(. Recently, Lasky and Donova examined the 



X 



stability and quasi- normal modes of black hole in TeVeS [15|. 

Additionally, it has been suggested how to distinguish TeVeS from GR observationally. For example, one could 
reveal the gravitational theory in the strong-field regime with the redshift of the atomic spectral lines emana ting from 
the surface of neutron star [lj], with the Shapiro delays of gravitational waves and photons or neutrinos 16], and 
with the spectrum of gravitational waves emitted from the compact objects (l7j . where the compact objects were 
assumed to be spherically symmetry. On the other hand, in this article, we will focus on the toroidal oscillations in 
the slowly rotating compact objects constructed in [l8j], which are associated with the r modes gravitational waves. 
The r modes arise due to the rotational effects and degenerate into zero frequency in the limit of non-rotatioin. 

The observations of stellar oscillations via gravitational waves are considered to provide a unique tool to estimate 
the stellar parameters such as mass, radius, rotation rate, magnetic fields, and equation of state (e.g., |19l - [23j ). which 
is called "gravitational wave asteroseismology" . The detailed analysis of the gravitational waves also makes it possible 
to determine the radius of accretion disk around supermassive black hole [24j or to know the magnetic effect during 
the stellar collapse [25| . 

In this article, as a first step to see the dependence of toroidal oscillations on the gravitational theory, we assume 
the Cowling approximation, i.e., we see only fluid oscillations and the perturbations of the other fields will be omitted. 
The more detailed study including the oscillations of the other fields will be done near future. This article is organized 
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as follows. In the next section, we review the fundamental parts of TeVeS and the stellar model with slow rotation 
in TeVeS. In Sec. IHIl we derive the perturbation equations describing the toroidal oscillations with the Cowling 
approximation and present the frequencies of toroidal oscillations as varying the stellar parameters. Finally, we make 
a conclusion in Sec. IIV1 In this article, we adopt the unit of c = G = 1, where c and G denote the speed of light and 
the gravitational constant, respectively, and the metric signature is (— , +, +, +). 



II. STELLAR MODELS IN TEVES 



A. TeVeS 

The details of TeVeS can be found in 16}. TeVeS is based on three dynamical gravitational fields; an Einstein metric 
g^, a timelike four-vector field , and a scalar field ip. Additionally there is a nondynamical scalar field a. The 
vector field satisfy the normalization condition with Einstein metric as g^U^U" = — 1 and the physical metric g^ is 
defined with the Einstein metric as 

9nu = e~ 2v g^ - 7UJA V sinh(2^). (2.1) 

The total action of TeVeS, S, contains contributions from the three dynamical fields mentioned the above and a 
matter contribution (see [6| for the details). This total action has two positive dimensionless parameters k and K, 
which are corresponding to the coupling parameters for the scalar and vector fields, respectively. Varying the total 
action S with respect to g^ v , U^, and tp, one can get the field equations for the tensor, vector, and scalar fields (see 
Q for the explicit field equations). Since the previous study about the neutron star structure in TeVeS has shown 
that the stellar properties are almost independent from the scalar coupling k [14j |. in this article we focus only on 
the dependence of vector coupling K. Although, the restrictions on K have not been discussed in great detail in the 



literature, in [14| they showed that to construct the stellar models K has to be less than 2 and also that K should 
be less than 1 to produce a realistic stellar mass. Thus, in this article we examine as varying K in the range of 
< K < 1. 



B. Slowly Rotating Relativistic Stellar Models 

As a background stellar model, we consider a slowly rotating relativistic star with a uniform angular velocity CI, 
where the rotational axis is set to be 9 = 0. Since the details for constructing such stellar models have shown in 
[l8[, in this section we describe only essential points. In the framework for slow rotation, we assume to keep only the 
linear effects in the angular velocity. Then, the stellar models are still spherical, because the deformation due to the 



rotation is of the order CI 2 . Those stellar models in TeVeS can be constructed by using the recipe shown in 14j and 
the metric in physical frame is given by 

d~s 2 = - e u+2 ^dt 2 + e^ 2ip dr 2 + r 2 e~ 2lf> (d6 2 + sin 2 6W0 2 ) - 2wr 2 e~ 2v sin 2 6dtd(j>, (2.2) 

where v, C, and lu are functions of the radial coordinate r. Up to first order of CI, the background fluid four-velocity 
of the star is described as 



(2.3) 



In TeVeS, another variable also needs to determine the rotational dragging, i.e., induced vector field due to the rotation, 
which is described as SU^ — —V(r) [l8|. With appropriate boundary conditions at the stellar center and infinity, 
by calculating two second order differential equations with respect to cj and V as shown in [l8|, one can determine 
the distribution of rotational frame dragging. About the stellar matter, we assume the perfect fluid described by the 
energy-momentum tensor 

f » v = [p + P) v?u v + Pg» u , (2.4) 

where p and P are the energy density and pressure in the physical frame, respectively. The adopted equations of 
state to construct the stellar model are the same ones as in [3j , which are polytropic ones derived by fitting functions 
to tabulated data of realistic equations of state known as EOS A and EOS II. The maximum masses of neutron stars 
with these equations of state in GR are M = 1.65M© for EOS A and M = 1.95M© for EOS II. That is, EOS A and 
EOS II are considered as soft and intermediate equations of state, respectively. 
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FIG. 1: Distribution of cj(r) in physical frame with different values of K for the stellar model with EOS A and Cl = 1 kHz, 
where stellar mass is fixed to be Madm = IAMq. In order to compare the results in TeVeS, the distribution in GR is also 
plotted with a solid line. 




As an example, the distribution of w(r) in physical frame for the stellar model with Madm = 1.4M© and f2 = 1 kHz 
is plotted in Fig. [TJ where Madm denotes the total Arnowitt-Deser-Misner (ADM) mass (see [3 for the definition 
of ADM mass). In the previous article, it was shown that the distribution of ui(r) depends on the value of vector 
coupling, K, as well as the total angular momentum of the star J, which is determined from the asymptotic behavior 
of uj(r) as uj(r) = 2.J /r 2 + 0(l/r 4 ) However, in this article we find that, if one would see the distribution of ui 
as a function of r / R where R is stellar radius, those distributions with different value of K are very similar to that in 
GR. 



III. TOROIDAL OSCILLATIONS 



In this article we focus on the toroidal oscillations with the Cowling approximation. That is, we consider only the 
fluid perturbation with axial parity and the other perturbations of scalar, vector, and tensor fields are neglected. We 
should notice that the Cowling approximation in GR is typically quite good for toroidal oscillations (axial parity) 
because this type of oscillations dose not involve the variation of density, while for spheroidal oscillations (polar parity) 
one can qualitatively discuss those frequencies but the error for typical relativistic stellar models could become less 



than 20% for the fundamental modes and around 10% for the pressure modes [26j. Similarly GR, we might expect 
the validity of the Cowling approximation for toroidal oscillations in TeVeS. However, as pointed out in |18| . for the 
slowly rotating star in TeVeS there exists the induced vector field due to the rotation and this induced vector field 
might play an important role in the toroidal oscillations, when one considers the metric perturbation. As a future 
work, we will examine the toroidal oscillations with the perturbations of the other fields. 

With the approximation of slowly rotating star as mentioned in Sec. UH we consider the rotational effect up to 
first order of fi. On this slowly rotating star, we add the fluid oscillations. Then, the leading order in the linearized 
field equations becomes ~ O(fie), where e expresses the order of fluid oscillations. Additionally, one can see that 
the coupling between the toroidal and spheroidal oscillation becomes higher order effects than 0(Cle) (see [28| for 
the ordering in GR). In this article, as a first step, we take into account the leading order of the fluid oscillations, 
in which we omit the coupling between the toroidal and spheroidal oscillations. Considering the toroidal oscillations, 
the Lagrangian displacement vector for the fluid perturbation is 

1 „ 1 



e = (?, e , = ( o, -z—a+Ytn, z—d 9 Y em ) , (3.i) 

V / V sine/ sine/ / 



where Z is a function of t and r, while Yg m — Yf m (0, 4>) is the spherical harmonic function. Then, the perturbations 
of four- velocity in the physical frame, Su^, can be written as 

5u r = 0, (3.2) 

5u 6 = -e- v - v ' 2 d t Z^—d^Y lmi (3.3) 
sins' 

5u* = e -^/ 2 d t Z^—dgY lm , (3.4) 
sine/ 

Su* = (Q-u) r 2 e-^- v sin 2 Su*. (3.5) 
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It should be noticed that in order to get the expression of Su*, we use the relation that g^u^Su 11 = 0, which is 
obtained the normalization condition for fluid four- velocity, i.e., g^, v u^u u — — 1. With the Cowling approximation, 
the axial perturbation of energy-momentum tensor is given by 

8 fu,u = fp + p^ (sfpff* + u»8u v ) . (3.6) 

The perturbation equation describing the toroidal oscillations can be obtained by taking a variation of the energy- 
momuntum conservation law, V^T^ = 0, which can be reduced as V„<5T'"' / — with the Cowling approximation. 
The explicit forms with fj, = 6 and <fi are 

a£ m {t,r)^—d (f ,Y em - fa m (t,r) cos0d e Y im = 0, (3.7) 
sm 6 

cos 

ae m (t,r)d e Y tm + Pi m (t, r)-— -dd,Y em = 0, (3.8) 

smt) 

where the coefficients ai m and /?£ m are 

a tm = dtZ + imfldtZ, (3.9) 

(3e m = 2 (w - fi) d t Z. (3.10) 



Calculating (Eq. (|3.7[) )/ sin6d^Y^ m + (Ea. p.8|l )9flY' e * n and integrating over the solid angle, one can get the single 
perturbation equation for the toroidal oscillations, such as 



dfZ + im 



n - — Q 



dtZ = 0, (3.11) 



where Y£ denotes the complex conjugate of Y( m and to derive Eq. (13. lip we use the following relations; 



2tt fir / -y 



deY; m d e Y e , m , + —^-d^Y t * m d^Y ( , m , sin OdBty = £(£+!) 5 mrn ,6 u , , (3.12) 
sm 



I I \ -r^ d eYi m d<t>Yt'm' ~ ^^rd^Y^dgYtm' ) sin 6 d8 d<j) = im6 mm , 5 u >, (3.13) 
Jo Jo \smf sinW / 

where Su> denotes the Kronecker delta. Now, assuming the time dependence of perturbation variable Z as Z(t, r) — 
Z(r)e lat > the frequencies of toroidal oscillations are given by 



£(£ + 1) 



(3.14) 



In spite of the fact that the adopted gravitational theory is different from GR, these equations (|3.11l) and (|3.14[) are 
same forms as in GR. That is, at least with the Cowling approximation, the background scalar and vector fields do 
not affect directly on the toroidal oscillations. 

In the Newtonian limit, in which uj — » 0, the spectrum of toroidal oscillations is discrete and the frequency becomes 
single value as a = [2/£/(£+ 1) — l]mO. However, when one considers the relativistic effect, i.e., frame dragging, the 
frequency of toroidal oscillations is not single value. Unlike Newtonian case, the frequencies become the function of 
u, i.e., the function of r. Namely, similar to GR case [28|, |29j, the frequencies in TeVeS could be continuous spectrum 
limited to a certain range. The allowed frequencies are determined with the value of ui inside the star. Since the value 
of co(r) is monotonically decreasing as r is increasing, the minimum and maximum values of ui inside the star can be 
determined those values at the stellar surface and center, respectively, i.e., uj(R) < u> < cj(0). Combing this evidence 
with Eq. (|3.14p . the frequency of toroidal oscillation should be limited in the range of 

|cr mi „| < |er| < |er 

max 

|, (3.15) 

where <7 m i n and a ma x are 

o-min = -mh + 2m (fi - u(R)) , (3.16) 



cr, 



max 



= — mQ 



1) 
2m 

£{£ + !) 



(fl-cj(O)). (3.17) 
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From observational point of view, it is not sure whether this continuous spectrum could be observed or not. In fact, if 
one considers in GR the fast rotating star (or the coupling with the spheroidal oscillations) , the frequencies of toroidal 
oscillations also show the discrete spectrum [13, So, the observed frequencies might be a single value in the range 
of Eq. (|3.15p . However, through the analysis as in this article, we believe that one can estimate how the gravitational 
theory affects on the frequency of toroidal oscillations. 




Q[Hz] 



FIG. 2: Frequencies with I = 2 in TeVeS with K = 0.45 are plotted as a function of the stellar rotational frequency Q for 
the stellar model with Madm = 1.4M© and EOS A, where the solid lines denote the values of a max , while the broken lines 
correspond to those of <7 m m. 

It is well-known that with the slowly rotating approximation, the frequencies of toroidal oscillations in GR are 
proportional to the angular velocity. On the other hand, the frequencies of toroidal oscillations in TeVeS with £ = 2 
are shown in Fig. [2 where the stellar mass for EOS A is fixed to be Madm = 1.4M Q with K = 0.45. From this 
figure, one can see that the frequencies even in TeVeS are proportional to il as same as in GR. Thus, we can define 
the new parameters, a max and a m i n , as 



and 



mQ, 



(3.18) 



which are independent of the value of O and to. Figs. [3] and [4] show these parameters in GR and in TeVeS as a 
function of I defined by 

H (3 - i9 » 

where I is constant and corresponding to the relativistic generalization of momentum of inertia for slowly rotating 
system [lH, |27|- In these figures, as increasing the stellar mass Madm along each line, the values of a max and a m ; n 
are also increasing. It should be noticed that the value of I has not observed yet, but I is an important parameter 
expressing the stellar configuration and that value might be determined if the phenomena of precession of rotating 
compact stars would be observed. From Figs. [3] and HI one can see that a max and a m i n are not so sensitive on the 
gravitational theory, while / depends strongly on the gravitational theory. In practice, if the values of a max and/or 
Omin would be determined, the value of I could change around 80% depending on the value of vector coupling K. 
Additionally, it could be found that the values of a max and a min are almost independent from the adopted equation 
of sate, but the corresponding value of I depends strongly on the adopted equation of state. That is, via the detailed 
observations of a max (or a m ; n ) and I, one could distinguish not only the gravitational theory in the strong-field regime 
but also the equation of state constructed the compact object. 

On the other hand, in Figs. [5] and El the ratios of a max and a m i n in TeVeS to those in GR, which are defined by 

(TcVcS) (TeVeS) 

r - Qmax fln H R ■ - Qmin C\1(W 

fmax — (GR) o-nu ymm — (GR) ' 
Omax a min 

are plotted as a function of the vector coupling K, where the stellar masses are fixed as Madm = 1.3Mq, 1.4M Q , and 
1.5M for EOS A, and Madm = 1.4M©, 1.6M Q , and 1.8M for EOS II. These figures show that with smaller K, 
the values of /3 m ax and /3 m i n become around 1 independently of stellar mass, while with larger K, one can see small 
deviation of a max and a m i n in TeVeS from those in GR at most a few %. Furthermore, it is also found that as the 
stellar mass becomes larger, the values of /3 max and /3 min are larger. So, with massive stellar model and with larger 
value of K, it might be possible to distinguish the gravitational theory by observing the detailed toroidal oscillations 
with the help of the observation of stellar mass. 
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FIG. 3: Values of a max and a m i n in GR and in TeVeS for EOS A as functions of I, where the solid lines correspond the results 
in GR and the dotted lines are those in TeVeS with different values of K, such as K — 0.2, 0.5, 0.7, and 1.0. 




FIG. 4: Similar to Fig. d but for EOS II. 



IV. CONCLUSION 



In order to examine the dependence of frequencies of toroidal oscillations on the gravitational theory, we consider 
such oscillations on the slowly rotating relativistic star in Tensor- Vector-Scalar (TeVeS) theory, where the angular 
velocity f2 is assumed constant. For this aim, we have derived the perturbation equations describing the toroidal 
oscillations of neutron stars in TeVeS and examine their specific frequencies as varying the vector coupling K and 
stellar mass, where as a first step we focus only on the fluid oscillations and omit the perturbations of tensor, vector, 
and scalar fields (the Cowling approximation). We find that the perturbation equations in TeVeS can be written 
same form as in general relativity (GR). Therefore, similarly GR, the frequencies of toroidal oscillations show the 
continuous spectrum at least in the frame with the Cowling approximation. Comparing the results in GR and in 



1.020 




K 



FIG. 5: Dependences of /3 max and /9 m in on K for EOS A, where the stellar masses are fixed to be 1.3Mq (traiangles) , IAMq 
(squares), and 1.5M© (circles). The solid and dotted lines correspond to the values of /3 m ax and fimin, respectively. 
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FIG. 6: Similar to Fig. [5] but for EOS II, where the stellar masses are fixed to be 1AM@ (traiangles), 1.6M (squares), and 
1.8M Q (circles). 



TeVeS, one can see that although it seems that those frequencies are not so sensitive on the gravitational theory, the 
value of relativistic generalization of momentum of inertia depends strongly on the gravitational theory as well as 
the adopted equation of state. So, observing the frequencies of toroidal oscillations and momentum inertia with high 
accuracy might be able to reveal not only the gravitational theory in the strong-field regime but also the equation of 
state in the higher density region. 

In this article, for simplicity, we adopt the Cowling approximation. That is, our consideration restricts only on the 
stellar oscillations. This means that we should do a more detailed study including the perturbations of other fields. 
With the oscillations of other fields, it could be possible to obtain the additional information about the different type 
of oscillation modes. That is, the observations of other oscillations can provide more accurate constraints on the 
gravitational theory in the strong-field regime. On the other hand, we should take into account the mode coupling 
between the toroidal and spheroidal oscillations, although such coupling is higher order effect. Considering this type 
of coupling, it is known in GR that the spectrum of toroidal oscillations could become discrete [3(|. The same can 
be expected in the case of TeVeS. If so, since the spheroidal oscillations depends strongly on the gravitational theory 



17j , via the mode coupling, the frequencies of toroidal oscillations could also depend strongly on the gravitational 



theory. Thus, the observations of toroidal oscillations might become more important to distinguish the gravitational 
theory. Furthermore, it might be important to study the dependence of magnetic effects on the toroidal oscillations. 
In practice, the quasi-periodic oscillations have observed during the decaying tail of giant flares and these phenomena 
are believed to be related to the oscillations of strong magnetized neutron stars [3l[ . Taking into account the magnetic 
effects, one might be possible to obtain the further constraint in the gravitational theory. 
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